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A method of using X-ray absorption spectroscopy (XAS) together with resolved grazing incidence geometry 
for depth profiling atomic, electronic, chemical or magnetic local structures in thin films is presented. The quan- 
titative deconvolution of thickness-dependent spectral features is performed by fully considering both scattering 
and absorption formalisms. Surface oxidation and local structural depth profiles in nanometric FePt films are 
determined, exemplifying the application of the method. 



I. INTRODUCTION 

Magnetic thin films have attracted a lot of attention due 
to their extremely high-density magnetic recording applica- 
tions 1 1 1. In this regard, a clear understanding of the macro- 
scopic magnetic properties requires a substantial knowledge 
of its dependence with layers thicknesses |2| and the com- 
plex microstructural effects frequently localized at the inter- 
face with the substrate or the surface of the films. Such effects 
can be studied using experimental techniques able to peer se- 
lectively in the depth of the films. In a previous letter 121 
we presented qualitative results using x-ray absorption spec- 
troscopy (XAS) with resolved grazing incidence (GI) to clar- 
ify the thickness-dependent magnetic properties in nanometric 
CoPt films. A depth dependent chemical order was revealed 
and the magnetic behavior was interpreted within this frame- 
work. In the present paper we provide a rigorous quantitative 
method for the deconvolution of the local atomic, chemical 
and magnetic structural depth profiles. This method is then 
illustrated by its application to oxidized FePt thin films @|. 
The proposed approach makes GI-XAS a unique tool to ad- 
dress the depth dependence of the local structural parameters, 
suitable for nanometric structures where this dependence is 
a crucial issue. Moreover this method provides a new venue 
to rigorously determine depth dependent electronic structure 
profiles using XANES (x-ray absorption near edge structure), 
which turns out to be crucial in understanding striking artifical 
interface materials |l5]|6l. 

Although the general phenomena of scattering and absorp- 
tion of x-rays by the condensed matter are nowadays quite 
well understood fT), they still are normally explored from 
unconnected viewpoints. A few well established techniques 
surpass this general rule with interconnected scattering and 
absorption techniques, providing invaluable additional selec- 
tivity compared to each approach used separately. DAFS 
(diffraction anomalous fine structure) gives site selectivity and 
local structural information IH. XAFS (x-ray absorption fine 
structure) extracted from reflectivity data gives local structural 
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information from surfaces and interfaces |9l. XSW (x-ray 
standing wave) locates impurities in bulk crystals and nanos- 
tructures using an interference field that provides spatial de- 
pendence to the x-ray spectroscopic yields from atoms within 
the field ifTOl . Similarly, glacing-incidence x-ray fluorescence 
(GIXRF) is a sensitive probe of chemical composition as a 
function of depth |11|. These techniques are based on sim- 
ilar approaches to the one presented here, however they are 
limited to near-perfect crystal structures, require well defined 
geometries and/or give a limited set of information. 



II. GRAZING INCIDENCE X-RAY ABSORPTION 
SPECTROSCOPY 

X-ray absorption spectra contain information about the 
ground state of the selected element in a material (local sym- 
metry, oxidation and spin states, spin-orbit coupling in the 2p 
and 3d orbitals, crystal field, covalence and charge transfer). 
As a matter of fact, in the case of 3d transition metals es- 
sentially structural information is obtained from the K edges, 
while more magnetic and electronic information is usually de- 
duced from L2.3 edge. XAS is not a surface technique by it- 
self, since the atenuation length of hard x-rays is of a few mi- 
crometers in any material. However, in the grazing incidence 
geometry near the critical angle for total reflexion, the x-ray 
beam is confined within a few nanometers from the surface. 
For this film studies, this confinement has the considerable 
advantage of minimizing the substrate contribution. 

The grazing incidence x-ray absorption measurements were 
performed at the Brazilian Synchrotron Light Laboratory 
(LNLS - Laboratorio Nacional de Luz Smcrotron). The setup 
includes 20 /zm-vertical slits limiting the beam size on the 
sample mounted on a high precision goniometer. XANES 
spectra were collected in the fluorescence mode at the D04B- 
XAFS 1 beamline 1121] with a Si ( 1 1 1 ) channel-cut monochro- 
mator. The incident beam intensity was monitored using a first 
ion-chamber. The reflected beam and fluorescence emission 
were simultaneously collected using a second ion-chamber 
and a 15-elements Ge detector, respectively. The fluorescence 
emission and/or x-ray reflectivity curves were used to cali- 
brate and select with an accuracy of w 0.01" the working graz- 
ing angle corresponding to a chosen penetration depth profile. 
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For an accurate energy calibration, the transmission through 
an Iron metal reference foil was monitored using a third ion- 
chamber 

The collected absorption spectra measured by the fluores- 
cence yield is a mix of contributions coming from different 
depths. To get quantitative information we must deconvolve 
them into their absorption contributions from each depth (z) 
into the films at each photon energy (E) and grazing angle (0). 
The electromagnetic radiation amplitude at each set of (E, 0, 
z) must be known to weightly sum the absorption contribu- 
tions as function of energy and angle {jiexpiE, 9)), as follows: 
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M,,p(£,0) = - / I{E,z,d)li{E,z)dz 



(1) 



where I{E,z,9) is the radiation intensity as function of E, z 
and 0; jj.{E,z) is the absorption spectrum contribution at the 
depth z and F is the normalizing factor I{E,z, 0)dz- The 
main difficulty to determine each iJ.{E,z) by solving this equa- 
tion is the initial calculation of I{E,z,9), which depends on 
how the layers structure of the film dynamically refract and re- 
flect the incident radiation, as a function of energy and depth. 
The formalism adopted to determine this intensity and the way 
to extract the depth dependence from GI-XAS spectra are de- 
scribed in the following section. 



A. Refracted and reflected amplitudes as a function of the 
penetration depth, photon energy and incident angle 

Several approaches EllIIllIIllIlIBlIIlllIimS 19 20 21 1 

can be used to estimate the transmissivity of x-rays inside a 
material. Those based on the dynamical diffraction theory 
give the most accurate results near the critical angle of to- 
tal external reflection. We deal here with conditions near the 
critical energy and angle for absorption and reflectivity reso- 
nances. Hence the method must include all dynamical reflec- 
tions and refractions conditions to determine the internal elec- 
tromagnetic wave amplitude in the samples. To fulfill these 
requirements we apply an approach analogous to the recur- 
sive Parrat's reflectivity method 121] to calculate the refracted 
and reflected amplitudes at every depth within a film formed 
by n layers, each one with different chemical contributions. 

Following the definitions by Parrat |21|, we consider an 
electromagnetic wave propagating into a material: 

E = £oe'["(^-^'-™l = g^^,l{i-5)Ck-7)-o>t]^-l3Ck-7) (2) 

where 5 and /3 are the real and imaginary part of the complex 
refraction index llTl fT4l l22l » = 1 - g - . 

The continuity at each interface between the media n and 
n — 1 of a film with layers gives the twin equation : 



a„-iE„^i +fl^^j£',f_j 



(3) 



{a„-iE„_i+a^\E^_i) f„_lki = {a„^E„ + a„E^) f„ki 

where /„ = {0^ ~ 25„ — 2ifi„yl^ for each n media and a„ = 
e ^ ^" " is a phase factor taking into account the absorption. 
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Figure 1 : (color online) Experimental refraction index h= 1 — 5 — / j3 
corresponding to FePt and Fe203. In (a) is the imaginary j3 compo- 
nent, and the real 8 component obtained with Kramers-Kronig trans- 
form is in (b). 



where d„ is the half penetration in the media n. E„ and Ej^ are 
the total and reflected electric field amplitudes in the media 
n. The solution of the equation [3] for the reflected amplitude 
Rn-i,n can be recursively determined by 
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The reflectivity at the interface between the air (or vac- 
uum) and the film is obtained after previous determination of 
Rn-i.ii at all others interfaces inside the film, considering that 
Rn-i.n = 0, as the infinitely thick substrate does not add any 
reflection. 

The electromagnetic radiation amplitude at each depth (z) 
inside a thin film can be determined in the same way by solv- 
ing the equations |3] to inform the amplitude £„ at each layer 
n and depth z- Isolating E„ in equation [3] and using the value 
Rn-\,n — '^n-i(^ff-i/^n-i)' '^^^ Straightforwardly obtains: 
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E„ can be recursively determined from the previous knowl- 
edge of the elements R„.„+\ and F„^\,„ calculated for the total 
reflectivity. Consequently, the amplitude at the upper interface 
of layer n is : 



bn^-^il+Rn,,+ l) 
GnEo 



(6) 



with Eq incident amplitude on the film surface (n = 0). 

The amplitude inside the layer « at an arbitrary position 

Zn relative to the top of the layer is then B„ = b„- e ^ ■'" " 
and the radiation intensity inside this layer n is I„{zi,,0,E) — 



b„ ■ 



.1 n^'H 



The total intensity I{z,9,E) at each depth z 

is given by the set of I„{zn,0,E) considering each thicknesses 
d„ and all possible n. The angle and energy dependences con- 
tained in the f„{8,l5,9) complex terms are fully mathemati- 
caly and computationally considered, where )3 and 5 are the 
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components of the refraction index, as shown in figure [T] for 
FePt and Fe203 compounds. 



B. Depth dependence of XAS spectra 

The electromagnetic radiation intensity I{E,z, 0) described 
above is used to determine the depth profile of XANES exper- 
imentally obtained in the grazing incidence geometry. This is 
performed by fitting XANES spectra for several grazing an- 
gles around the critical angle, considering the x-ray attenua- 
tion inside the material. The result of this process is the strati- 
fication in layers (of thickness dz) of the XANES information. 
The XANES spectra for each depth are fitted as a linear com- 
bination of reference spectral contributions previously deter- 
mined. As the structure of layers inside the film can change 
dynamically in the fitting process, the I{z,0,E) intensity must 
be calculated at each self-consistent fitting iteration. 

To determine the depth dependence of jiexp, it must be 
found a set of iJ.{E,z) data that when convoluted with 
I{E,z, 9) in the equation[T|simultaneously fits Hexp{E, 9) mea- 
sured for several 9. This is more easily done rewriting the 
equation [T] in a discrete form considering the sum in z with 
steps of Az in depth: 

M«p(£,0) = ^f;/(£,Zi,0)M(£,zO (7) 

1=0 

where zo = 0, Zoo = °o, Zi+i -Zi = Az and T = YT=QKE,Zi, 9). 

It is easily seen that the proportional contribution (PC) for 
each Az layer at depth z, for the signal pLexp{E,9) is deter- 
mined by PC{zi) = I{E,Zi, 9)/r. If the experimental spectra 
can be considered as a linear combination of several indepen- 
dent contributions of XANES features, jigxpiE, 9) can be con- 
sidered as a linear combination of q spectral contributions (j) 
each one weighed by a factor w,-: 

fiexp{E,9) = ^f^fij{E)-Wj{9) (8) 

where 11 is the normalization factor 11 = Ly=i 

Considering the equations [t] and |8] each n{E,Zi) can be 
written as pL{E,Zi) = n ■Wj{zi) where n(z;) = 

'L]=iWj{zi). Therefore: 

(=0 7=1 

The objective of determining all absorption contributions 
{lXj{E)) and its equivalent weight {wj{zi)) for each depth z, 
can be achieved by fitting the experimental spectra pLexp{E, 9) 
with equation |9] Is important to note that since ii {E) for 
each spectral contribution at each layer is included on both 
I{E,z,9) and ix{E,z) of equation |7j the deconvolution of 
IJ.exp{E, 9) must be a self-consistent procedure in terms of j3. 
In other words, I{E,z, 9) must be computed at every iteration 
of the iJ.exp{E, 9) fitting procedure. 



III. APPLICATION: DEPTH PROFILE IN A FEPT 
MAGNETIC THIN FILM 

The FePt film studied here was grown by sputtering from 
pure targets of elemental Fe and Pt. It was deposited on 
MgO(lOO) substrate with a pre-deposited Pt fcc(lOO) buffer 
layer The substrate temperature was kept at 500°C to ensure 
a good chemical order and perpendicular magnetic anisotropy 
im. The composition and thickness of the film was checked 
by Rutherford Backscattering Spectroscopy (RBS) confirm- 
ing the equiatomic ratio (51% Fe and 49% Pt) and a thickness 
of 103 nm for the FePt layer. 

When deposited at high substrate temperature (> 400°C) 
FePt thin films without cap layer protection are easily oxidized 
materials 11231 . The presence of a surface oxidation is clearly 
observed in the GI-XAS measurements at the smallest graz- 
ing angle presented on figure [2] and the oxide is identified as 
Fe203. The contribution of the oxide layer decreases rapidly 
for increasing angles, indicating that this layer is limited to a 
few A. This is clearly seen by the decreasing (increasing) fea- 
ture at 7133 (7115) eV and the shift of the spectra to lower 
energies presented in figure |2] (upper inset). 

Simple tabled or calculated 5 and /3 components of the re- 
fraction index n — 1 — 5 — ;j3 [ T4ll22]| cannot be used to cal- 
culate I{E,z,9), given that it would include approximations 
not vaUd when the spectral features near a critical energy (ab- 
sorption edge) are the desired information. The imaginary )3 
component of the experimental refraction index was obtained 
for reference FePt and Fe203 samples using their absorption 
spectra scaled to tabled absolute values far from the absorp- 
tion edges |fT4l l22l . as shown in figm-e[T^. Kramers-Kronig 
transforms 1 24 25l l26l l27l l28l l29l l30ll were used to determine 
the correspondent real 5 component (fig. [TJi). Figure [3] shows 
the simulated intensity I{E,z, 9) for a model thick FePt film 
with a flat 5 nm Fe203 layer on its surface, using equation |6] 
and the refraction index data for each layer presented in fig- 
ure [T] The XANES structures must definitively be taken into 
account when determining the experimental refraction index 
in order to include all energy/angle dependences in I{E,9), 
as is ratified by the strong non-linear dependence on I{E) for 
different grazing angles shown in fig. [3^. The intensity I{9) 
at E= 7130eV for representative z values is shown in fig. [3]3, 
which exemplifies multiple reflection interference effects near 
the critical angle for each penetration. The interference con- 
dition for the first Fe203 layer (5 nm) is drastically different 
from the resonances for the deeper layers, where the contribu- 
tions arise essentially from the FePt alloy. These simulations 
emphasizes the strong need to consider all dynamical reflec- 
tions at the interfaces to accurately calculate I{E,z, 9). 

On the other hand it is worth noting that although correc- 
tions due to fluorescence self-absorption effects might be im- 
portant in some cases, these effects are not significant for the 
angular range and penetration depth discussed here (<1% in 
the final error bars). 

Equation [9] considering the iteratively determined 
I{E,z,9), was used to simultaneously fit the experimental 
XANES spectra taken at several grazing angles (fig. [2]i. 
Reference XANES spectra of Fe203 and FePt were used in 
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Figure 2: (color online) Fitted XANES experimental data as a func- 
tion of the energy and several grazing angles for a FePt film. The 
best model to fit the experimental data consists of a gradient between 
a oxidated surface with pure Fe203 layer (0.2 nm thick), with inter- 
mediary layers of both oxide and metal (0.7 and 0.5 nm thick), to the 
pure FePt metal inside the films. Upper inset : GI-XAS measure- 
ments for increasing grazing angles showing the dependence of the 
XANES features as a function of the penetration depth, or grazing 
angle. The arrows correspond to an increase in the grazing angle. 
Lower inset: weight function w for the two components (FePt and 
Fe203), for the oxidation profile corresponding to the best fit. 
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Figure 3: (color online) Calculated intensity of 

the electromagnetic radiation inside a thin film of 
[Fe2O3]5„in/[FePt]y5„n,/[Pt]50nm/MgO. In (a) is the intensity 
as function of photon energy and grazing angle, at fixed penetration 
of 2 nm. In (b) is the dependence with incident angle and the depth 
penetration, with fixed energy of 7130 eV. 



the fit. 

The several experimental spectra simultaneously fitted, 
shown in figure |2] enable us to determine the complex layers 
structure in the depth profile, beyond a simply oxidized thick- 
ness determination. Different models of layers structure were 
considered for the oxidized FePt surface. Although a flat top 
oxide layer is well-suited to illustrate the general behavior, it 
is by far not the right solution to fit our data. The depth profile 
analysis shows clearly that the oxide not only covers the FePt 
film but penetrates beneath the film giving rise to a fraction- 
ated buried layer composed of the oxide and FePt. The best 
fit modelturns out to be a gradient between a thin oxidized 
surface with pure Fe203 layer and intermediate layers of both 
oxide and metal alloy down to 1 .4 nm from the surface. The 
top 100% oxide layer is 0.2 nm thick, followed by two mixed 
layers with 70% and 30% of Fe203 and 0.7 and 0.5 nm thick, 
respectively. The weight function w for each component is 
shown in the lower inset in figure |2] 

It has been reported |23| that for FePt films an Fe oxide 
layer would form on the surface due to Fe migration to the 
oxide/metal interface during the growth at high temperatures. 
As a result, there might be a composition variation with in- 
creasing film depth. In the metallic layer just below the ox- 
ide, Fe content should be lower than that of as-deposited film 
while Pt content should be higher. Our results confirm a com- 
positional variation over the film depth, but supports a more 
complex picture. As known from literature, sputtered metal- 
lic films have some tendency to form pillars. We interpret the 
gradient as resulting from the decoration of these pillars by the 
oxide that fills the empty space between pillars and oxidizes 
the very interfacial Fe atoms, rather than a continuous rough 
surface. 

We should finally include an additional remark about the 
resolution of the depth profiles. Due to the exponential de- 
cay of the radiation intensity inside the film, the depth probed 
and final resolution of the method are intrinsically correlated 
and strongly dependent to the contrast between the refraction 
index of each layer material in the film. For instance, in the 
Fe203/FePt case the refraction index of Fe203 is factor three 
smaller than for the FePt material. In this case the profile vari- 
ation is confined within 2 nm near the surface and the depth 
resolution is of order of one angstrom. If the compositional 
gradient were deeper into the film, the profile resolution would 
be lower for the internal layers. 



IV. CONCLUSION 

Scattering and absorption phenomena are intrinsically in- 
tercorrelated when grazing incidence reflection and refrac- 
tion are combined to x-ray absorption spectroscopy. The ap- 
proach of GI-XAS presented in this article fully considers 
both scattering and absorption formalism for depth profiling 
the atomic, electronic, chemical or magnetic local structures 
in thin films with nanometric resolution. This formalism, not 
facing intrinsic limitations or approximations, can be applied 
to deconvolve the depth dependencies of not only XANES in- 
formation as exemplified here, but also XRF (x-ray fluores- 
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cence) and XMCD (x-ray magnetic circular dichroism) sig- 
nals in the fluorescence or reflectivity channels from thin films 
and multilayers. 
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